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ABSTRACT 

Priedrichs,  in  his  studies  of  difference  approxi- 
mations to  solutions  of  symmetric  hyperbolic  systems, 
has  shown  that  difference  schemes  with  positive 
coefficients  are  stable.   In  this  paper  we  show  that 
this  stability  criterion  of  Priedrichs  is  valid  for 
difference  approximations  to  arbitrary  hyperbolic 
systems.   The  proof  relies  on  the  convexity  and  mono- 
tonic  ity  of  eigenvalues  as  matrix  fvmctions  over  a 
linear  space  of  matrices  with  only  real  eigenvalues. 
These  theorems,  well  known  in  the  symmetric  case,  are 
proved  with  the  aid  of  theorems  concerning  the 
dependence  of  solutions  of  hyperbolic  equations  on 
their  initial  data. 
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NYO-797i; 

DIFFERENTIAL  EQUATIONS,  DIFFERENCE 
EQUATIONS  AND  MATRIX  THEORy''^+ 

1»   Intro  due  ti on «, 

In  this  paper  we  shall  study  sets  of  real  matrices  all 
of  whose  linear  combinations  with  real  coefficients 
have  only  real  eigenvalues*   We  shall  show  that  over  a 
linear  space  spanned  by  such  matrices,  the  largest  eigen- 
value is  a  convex  function,  the  smallest  is  concave,  and 
all  eigenvalues  are  monotonic  with  respect  to  the  partial 
ordering  induced  by  the  convex  cone  of  matrices  with 
positive  eigenvalues.  This  is  a  generalization  of  classical 
results  known  for  symmetric  matrices  -  a  special  but 
important  subset  of  our  matrices.   Matrices  of  this  kind 
occur  as  coefficients  of  first  order  vectorial  hyperbolic 
equations.   The  proofs  of  the  matrix  theorems  mentioned 
before,  as  well  as  some  related  ones,  are  based  on  known 
theorems  in  the  theory  of  hyperbolic  equations  and  will 
be  presented  in  section  3«>   The  relevant  portions  of  the 
theory  of  hyperbolic  equations  are  sketched  in  section  2» 

*  This  work  was  carried  out  xinder  Contract  AT(  30-1  )-l480 
■  with  the  United  States  Atomic  Energy  Commission, 

-J"  Presented  to  the  American  Mathematical  Society  in 

Rochester,  December  1956  and  Washington,  October  1957« 
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In  section  I|.  we  apply  one  of  the  matrix  theorems  of  section 
3  to  derive  a  criterion  for  the  stability  of  difference 
schemes  used  to  approximate  solutions  of  differential 
equations.   This  criterion  is  an  extension  from  the 
symmetric  to  the  nonsymmetric  case  of  Priedrichs' 
stability  theorem  for  difference  equations  with  positive 
coefficients 6  Applications  to  various  specific  difference 
schemes  are  given  at  the  end  of  the  section. 

We  close  this  section  by  constructing  a  pair  of 
matrices  with  the  aforementioned  property  yet  which 
cannot  be  simultaneously  made  symmetric  by  a  similarity 
transformation.   Our  example  will  result  from  pertiirbing 
the  pair  of  symmetric  matrices 


^o  = 


10    0 


0-10 
\0   0    0^ 


Eveiry  linear  combination  aA  +  bB   with  a  and  b 
•^  o     o 

2    2 
real  and  a  +  b   =1  has  the  eigenvalues   -1,  0  and  1, 

Since  the  set  of  matrices  with  real  and  distinct  eigen- 
values is  open,  it  follows  that  if  X  is  any  (real) 
matrix,  then  for   e   small  enough  all  linear  combinations 
of  A   and  B,   vdiere   B  =  B  +  eX,   also  have  real  and 
distinct  eigenvalues.   We  shall  show  now  that  with  an 
appropriate  choice  for  X,   A   and  B  cannot  be  made 
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symmetric  by  the   same   similarity  -cransformationj.      Suppose 
to   the  contrary  that   there   exists   a  matrix     S      svich  that 
both     S"  A^S     and     s"  BS      are   ssnnrnetric.      Let     U     be   the 
unitary  matrix   i'.'hich  diaji'onalizes     S     A-S;    denotr^     SU     by 
T     and     T'-^'A^T     by     A.     Both     A        and     A      are   diagonal 
and,   bein;:^  similar,   have   the   same  eigenvalues,;      Therefore 

o  o 

Multiplying  by     T     we  get     A^T   =  TA^ ,      i.e.      T     ana      A^ 

coiTiniute*      Therefore     T     must  have   the   came  el'-:;envectors 

as     A^,      i».e*      T      is   diar-:onal-.'      On   the  other  hand 
o 

T-'^B?  =  u'-'-.'r^BSi'     is   also   SY/.aoetric .      Since     T     is 


4- 


dia,''^ona?(. ,    this   impiles,    it   Is.   easy    to   check »    tha 


^?1^ 32^13  ~  ^1?^'''^^31      * 

r.ixt  we  can  easily  cnoose  X  so  that  tnis  condition  is 
not  satisfied* 


1  Strictly  speaking ,  the  diagonal  elements  in  A  and 
A   may  not  appear  in  the  same  order,  but   A  can 

always  be  brougjat  into  A   by  a  unitary  transforina' 
tion* 
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2>   Hyperbolic  SquatlonSo 

We  shall  deal  with  first  order  systems  of  partial 
differential  equations  of  the  form 

m       , 

(2.1)  u^  ^fe'  V^'"'"^^ 

where  u  is  a  vector  of  n  unknown  functions,  and  A 
are  real  constant  coefficient  matrices.   The  initial 
value  problem  is  to  find  a  solution  of  this  equation 
with  prescribed  values  at   t  =  0  : 

(2.2)  u(x,0)  =  <^(x)   V 

The  initial  value  problem  is  said  to  be  properly  posed 

if  (2.1),  (2*2)  has  a  unique  solution  for  all  sufficiently 

--  say  k   times  —  dif ferentiable  data  ^(:-:}. 

We  start  with  the  following  classical  result: 

Theoren.  2»1 :   vixe  initial  value;  pL^oblem  ( 2s  1 ) ,  ( 2 ♦2  ) 
is  properly  posed  if  and  only  if  all  linear  combinations 
A  *  ?  of  the  coefficients  A   with  real  coefficients  £,, 
have  only  real  eigenvalues.   The  solution  depends  continu- 
ously on  the  data ,  i . e ^,  the  value  of  the  solution  at  any 
point  is  a  continuous  linear  functional  of  the  data  in 
the   C  -topology  for  N  large  enoughs 

An  equation  whose  coefficients  satisfy  this  condi- 
tion is  called  hyperbolic. 

The  necessity  of  this  condition  has  been  demonstrated 
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by  Hadamard  [8]  and  Banach  [1],  page  kk"!*  Their  arguments 
have  been  extended  by  the  present  author  to  equations  with 
variable  coefficients,  [16] » 

The  solution  can  be  constructed  either  by  Fourier 
synthesis  or  by  the  method  of  plane  waves  (see  Fritz  John, 
[13]  #  or  Cotirant-Lax,  [3])«  Since  the  proof  is  usually 
oai»ried  out  under  sllglitly  more  stringent  conditions,  we 
sketch  the  proof  very  briefly J 

Suppose  that  the  initial  function  ^     depends  only 
on  a  given  linear  combination  of  the  space  variables s 

<^(x)  =  ^(y),   where  y  is  the  scalar  quantity 
X  • €,  5  some  (real)  vector  in  the  dual  of  x-space»  We 
shall  construct  a  solution  u  of  (2*1)  with  initial 
value  \lf  which  Is  a  i'unciion  or  j     and   t  only»  Denote 
u(x,t)  by  v(y,t)s  we  have  the  I'ollowing  equation  for 

VI 


where  A  =  >   ^^  ^  «• 

To  solve  this  equation,  introduce  new  iniknowns 

'^      -1 
w  =  Tv,  T  so  chosen  that  A  =  TAT    is  lower  triangular, 

i,eo   a,  -  =  0  for  i  <  /£  «  According  to  known  theorems 

of  matrix  theory,  such  a  matrix  T  exists.  It  can  be 

chosen  to  be  real  if  the  eigenvalues  of  A  are  real. 
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and  chosen  so  that  the  magnitude  of  T  and  T"   (i»e,  that 
of  their  elements)  does  not  exceed  a  quantity  that  depends 
only  on  the  dimension  of  the  space.   E«g» ,   T   can  be 
chosen  to  be  an  orthogonal  matrix. 
The  equation  for  w  is 

(2.1*  )  w^  =  AWy  , 

Since  A  is  lower  triangular,  the  first  equation  contains 
w   only  and  has  a  unique  solution  once  w   Is  specified 

on  the   y   axis.   The  second  equation  contains  w   and 

2  1 

w   onlyj  substituting  the  already  knoxvn  value  of  w   we 

2 

can  determine  w   by  integration  along  the  second 

characteristic*   Recursively  we  can  determine  all 
components  and  can  give  a  boxind  for  the  magnitude  of  all 
components  and  their  first  derivatives  in  terms  of  bounds 
for  the  first  m+1  derivatives  of  the  data  \}f  independ- 
ently of  ^  » 

Given  an  arbitrary  but  smooth  function  (|i   of  x  , 
,..,  X^,   it  can  be  represented  as  a  superposition  of 
"plane  waves,"  i»e<.  of  functions   '!(x  £,)      (   e  John,  [13]  )• 
If  4   is  smooth,  the  functions   v]/   ai't;  ^a.icewlse  smootn, 
say  m+1  times  diff erentiable.   It  follows  then  from  our 
previous  estimate  that  the  superposition  of  the  correspond- 
ing solutions   v(x»4>t)   exists  and  has  continuous  first 
derivatives. 
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The  much  more  delicate  question  of  when  Cauchy's 
problem  is  correctly  posed  for  equations  not  homogeneous 
in  the  highest  order  term  is  discussed  in  Garding,  [?] 
and  A«  Lax,  [ll+lo 

Take  a  point   P  outside  of  the  initial  hyperplane* 
It  is  well  known  in  the  theory  of  wave  propagation  that 
the  value  of  u  at   P  is  uniquely  determined  by  the 
initial  data  on  a  compact  portion  of  the  initial  hyper- 
plane  (finite  signal  speed).  We  shall  give  now  a  precise 
description  of  this  phenomenon? 

The  value  of  u(P)  is  a  linear  functional  depending 

I  k 

continuously  on  the  initial  data  <J  in  some  C   topology* 

Therefore  u(P)  can  be  represented  as  a  scalar  product 

u(P)  =  ((t,R)   , 

where  R  is  some  distribution.   It  is  well  known  in  the 
theory  of  wave  propagation  that  the  support  of  R 
called  the  domain  of  dependence  of  P  --  is  compact* 

This  important  property  of  solutions  can  be  expressed 
without  any  reference  to  a  representation  as  follows: 

Theorem:   Given  any  point   P  there  exists  a  smallest 
closed  set  K  =  K(P)  on  the  initial  hyperplane   t  =  0 
such  that  if  the  initial  values  4(x)  of  a  solution  u 
of  fe»l/  vanish  in  any  open  set  0  containing  K^   then 
u(P)  =  0. 
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To  prove  this  theorem  we  have  to  show  that  if   |k| 
is  a  collection  of  closed  sets  with  the  above  property  — 
and  from  now  on  we  shall  call  these  sets  of  uniqueness  -- 
then  their  intersection  is  also  a  set  of  -uniqueness.   We 
do  this  in  two  steps, 

Leimna  2.1;  The  intersection  of  two  sets  of  unique- 
ness is  again  a  set  of  uniqueness. 

Proof;   We  have  to  show  that  every  solution  u  of 
(1.1),  whose  initial  value  u(x,0)  =  4»(x)  is  zero  in  an 
open  set  0   containing  K,  H  Kp,   vanishes  at  P.   Since 
the  closed  sets  K, /I  0   and  Kp/1  0   are  disjoint  we  can 
construct  two  open  sets   0,   and  0-  containing  K,  H  0 
and  Kp/IO   respectively,  so  that  their  closure   0,   and 
Op  are  also  disjoint.   Next  we  construct  an  infinitely 
differentiable  function  p(x)  which  is  equal  to   0  in 
0, ,   equal  to   1  in  Op   and  smooth  everywheroo   Define 
<t^  and  c^p  as 

<J^  =  p<i   ,     ^^  =   (l-p)<i   , 

and  denote  by  u,   and  Up  the  solutions  with  initial 

values  ^,      and  (^p  respectively.   The  existence  of  u, 

(P ) 
and  Up  follows^  '  from  Theorem  2.1,'  furthermore,  since 

solutions  are  uniquely  determined  by  their  initial  data> 

it  follows  that  the  solution  u  with  data  ^      is  equal 

2  At  this  point  we  use  the  hyperbolic  character  of  {2.1}. 

-  11  - 


to  u,  +  U2  • 

Clearly,   (^,    vanishes  in  0,  U  0  and  ^p  ^^ 
Op f  0  ;   since  each  of  these  open  sets  contains  a  set  of 
iiniqueness   (K,   and  Kp  respectively),  it  follows  that 
u,(P)  and  U2(P)   are  zero;  but  then  so  is  u(P)  = 
n^{T)   +  U2(P). 

Repeated  application  of  Lerama  2.1  shows  that  like- 
wise the  intersection  of  a  finite  nxnnber  of  sets  of 
uniqueness  is  again  a  set  of  imiquoness. 

Lemma  2«2:  The  intersection  K  =  /IK.  of  any 
number  of  sets  of  uniqueness  K .  of  which  at  least  one 
is  compact  is  again  a  set  of  Tiniqueness» 

Proof;   Let  0  be  an  open  set  containing  K  =  fl^^" 
By  the  Heine-Borel  Theorem,   0   also  contains  the  inter- 
section of  a  finite  nimber  of  the  sets  K.  •   Since  a 
finite  intersection  was  already  shown  to  be  a  set  of 
viniqueness,  the  conclusion  follows. 

Thus  we  have 

Theorem  2.2;  Every  point  P  has  a  smallest  set  of 
xmlqueness  K  on  the  initial  hyper plane   t  =  Oj   this 
set  is  called  the  domain  of  dependence  of  the  point  P  • 

Prom  now  on  we  take  the  point  P  to  be   (0,1), 
where  0  denotes  the  origin  on  the  initial  hyperplane. 
Its  domain  of  dependence  K  is  a  pointset  uniquely 
determined  by  the  equation  under  consideration, 
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Whereas  the  precise  determination  of  this  pointset  is  a 
delicate  matter,  it  is  easy  to  determine  its  convex  hull 
(i.e.,  the  smallest  convex  set  containing  it)  in  terms  of 
the  coefficients  A   of  equation  (2.1  )«► 

We  recall  the  notion  of  the  support  ftinction  h(^) 
of  the  pointset  K  ; 

(2.3)  UV   =  Sup  x€ 

xeK 

Here  §  is  any  (real)  point  of  the  dual  of  x-space 
and   x5   denotes  the  scalar  product  >  ^^a   • 

As  is  well  known,  the  support  function  of  a  pointset 
is  uniquely  determined  by  its  convex  hull  and  vice  versa. 

Theorem  2.3»   The  support  function  of  the  domain  of 
dependence  K  of  the  point   P  :   (0,1)  in  equation  (2.1) 
is 

(2.1+)  ^(^)=^naxA*^ 

where     X^      (X)      denotes   the   largest   eigenvalue  of   a 
matrix     X     and     A  *  ^     abbreviates     ^       A-'^ .    , 

Proof;     We   shall  prove    {2»k)   for      ^2.  ~  '^*    ^2  ~  ^y' 
£     =  0    •     We  have   to    show   that  for  any  positive      e      the 
half  space     x     <Xn      (A)-e     is  not   a  set  of  uniqueness 
whereas  the  halfspace     x     <  Vnax^^    )   "•"  e     is* 

To   show   the  first  half  of  this  proposition  we  have 
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to  exhibit  a  solution  whose  initial  values  vanish  for 
X  <  Ti    -  e  but  which  does  not  vanish  at  P  •  To 
construct  such  an  example  we  introduce  new  unknowns 
V  =  Tu  ,   T   so  chosen  that  A  =  TA'T?"   is  lower 
triangular  i.e.   a^^ /^  =  0  ^<^^     i  <  /.   For  v,   equation 


(2,1)  reads 


."t-i 


X 

Our  choice  for  v  shall  be  a  function  that  does  not 
depend  on  the  variables  x  ^...^x  .  So  the  equation 
for  V  reads 


v^  =  Av^^   , 

where  A  =  TA  T"  is  lower  triangular.  Componentwise, 
the  equation  reads 

(2.1')  vj=^aj^v^^  . 

The  diagonal  elements  of  A  are  the  eigenvalues 
H»**"^n  °^     ^       ^^   ^°^^   order;  say  /^^  is  the   j   ; 

^jj  =  ^max  • 

Choose  the  initial  values  of  all  but  the  jth  component 
of  V  to  be  zero  on  the  whole  initial  line,  and  choose 
v-'Cx,©)  to  be  different  from  zero  in  some  small  neigh- 
borhood of  X  =  \  ,  and  equal  to  zero  outside  this 
nei^borhood.      Solving   the  first      j-1      equations   in 
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(£•1  )  recursively  we  conclude  that  the  first   j-1 
components  of  v   are  zero  for  all  x    and   t^   There- 
fore  the   j    equation  reads 

which  asserts  that  v^  is  constant  along  every  line 

with  slope  -  X    *   In  particular,   v*"!  P)=v*'(0,l  )= 
■^     max 

v^ (7.^,0)  j^   0  . 

To  prove  that  the  half  plane  x  <  ^^^^^^  "•"  e  is  a 
set  of  uniqueness  we  rely  on  the 

Holitipyen  Uniqueness  Theorem;   Let   D  be  a  domain 
in  the  hyperplane   t  =  0   and  let   D(-f)j   be  a  smooth 
deformation  of  D  in  x,t   space  keeping  the  boundary 
points  fixed ^-^   Suppose  that  none  of  the  surfaces 
"DiV)      is  characteristic  at  any  point;  then  D  is  a 
set  of  uniqueness  for  every  point  swept  out  by  the 
family  of  surface Sa> 

We  recall  that  a  surface  is  characteristic  at  a 
point  if  the  direction  cosines  n^yJ^t  o^'  i'ts  normal 
satisfy  the  characteristic  relation 

det(n^I  -  ^ZZ  ^  k  ^^^  "^  ° 

or,  in  other  words,  if  n^  is  an  eigenvalue  of 

3  I.e»  a  continuously  differentiable  mapping  of  the 
cylinder  xeD,   01^<1  i^^o  x,t-space, 
(x,0)  going  into  x  = 'x  7  t  =  0  and  which  is 
independent  of  T  on  the  boundary. 
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A»  n  • 

To  apply  Holmgren' s  Theorem  to  the  situation  at 
hand,  choose  D  to  be  the  domain  which  contains  all 
points  cut  out  of  the  half -space  x  <  \         +  e  by  the 
sphere   |x|  =  R  • 

Take  the  case  that  X_    is  positive;  (if  not  the 
change  of  variables  x  =  x  -  at  will  accomplish  this). 
We  define  the  deformed  surfaces  DCtt)  ^J 

X 


t  =  r(i  -  ^-^)   (1  -  -IfJ-)  ,  0  <  r <  1  , 

ax 


For  R  large  enough,   D(r)  is  nowhere  characteristic 
for  0  <  "C  <  1  ♦   For  we  have 

n  .  =  0(g)  f  j  =  2,  •••,   m  • 

Clearly,  the  matrix 

I £— A^ 

"K       +e  * 
max 

is  nonsingular  for  e  positive  and  'T^    between  zero 

and  one,  by  definition  of  X     as  the  largest  eigen- 

*'  max  ° 

value  of  A  .  But  since  the  set  of  nonsingular  matrices 
is  open,  the  characteristic  matrix  is  also  nonsingular 

I  -  A  *  n 
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for  R  large  enough. 

Since  the  point   P  =  (0,1)   lies  on  the  surface 
D(l),   the  proof  of  theorem  2«3  is  complete. 

Another  proof  of  theorem  2,3  can  be  given  on  the 
basis  of  an  explicit  representation  for  the  value  of 
u  at  P  in  terms  of  the  initial  data*   Such  a  formula 
can  be  derived  as  folloi«:s: 

Define  the  "influence  function"  R  as  the  solution 
of  the  backward  initial  value  problem 

L''  R  =  0 

R(x,l)  =  6(x)  I   , 

where  L   denotes  the  adjoint  of  the  operator  L 
occurring  in  (2»l)o  This  initial  value  problem  can  be 
solved  explicitly:  R  is  the  Fourier  transform  of 

i(l-t)A*»§   . 
e 

Apply  Green's  forraula 

t=l 


/  /    (R«>Lu  -  u  *l'^R)    dx  dt     =   /    u  'R   dx 
0  gm  ^ 


t=0 


to  any  solution  u  of  (2«1)  and  to  the  influence 
function  R  ;  we  obtain 

u(P)  =  /<Kx)  R(x,0)  dx   • 

The  Po\arier  inverse   of     R(x,0)      is      e    "        ^    ',      since  we 
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have  assumed  that  all  eigenvalues  of  A  •  ^   --  aiid  thereby 
those  of  its  adjoint,   A  •  ?   --  are  real  for  ?  real, 
the  function  e"""   '  ^  grows  at  most  polynomially  for  ^ 

1  A*'''  •  P 

real.  For  complex  values  of  ^  ,   e     ^  is  entire 
analytic  of  exponential  growth.  Therefore,  according  to 

Plancherel  and  Polya' s  extension,  see  [I?],  of  the 

IA'^  •  E 
Paley-Wiener  theorem,  the  Poiirier  transfonn  of  e     ^ 

has  compact  support  K,   and  the  support  function  h(»7) 

of  the  pointset  K  is  equal  to  the  rate  of  growth  of 

A*  -^ 
e      •  This  brings  us  back  to  formula  (2»!4.)» 
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3«   Some  Matrix  Theorems, 

Theorem  3»ls   Let   {xj   be  a  linear  space  of  n  bj 

n  matrices  over  the  reals  such  that  every  X  has  only 

real  eigenvalues.   Then  \     ., ( X ) ,   the  largest  eigenvalue , 

is  a  convex  matrix  funotionc   Likewise   A  .  (X),   the 
mm    *   

smallest  eigenvalue.  Is  a  concave  matrix  function* 

Proof:   Let  A  ,  A  ,».,A   be  a  generator  of  the 
space   [x}  •   The  differential  equation  (2»1)  with  these 
matrices  as  coefficients,  is  hyperbolic  according  to  • 
Theorem  2«1»   Therefore,  the  concept  of  the  domain  of 
dependence  K(P)  ,   P  =  (0,1),   can  be  defined}  according 
to  Theorem  2.3,  the  support  function  h(^)   of  K(P)  , 
P  =  (0,1)   is  V^(X)  ,   X  =  YH   ^i-A*'  '-      But  the  support 
ftmction  of  any  domain,  b  being  defined  by  (2«3)  as 
the  supremum  of  linear  functions,  is  convexo   This 
completes  the  proof. 

In  the  symmetric  case,  this  result  follows  from  the 
classical  maximum  principle  for  the  largest  eigenvalue: 


^.a^"^>  =|||y|P=i  (^•''^' 


The  concavity  of  \.    (X)   follows  from  the  relation 

\.    (X)  =  -  X    (-X)  o   From  this  it  follows  that  the  set 
Dim        max 

of  positive  (non-negative)  matrices  X,   i.e.,  matrices 
with  positive  (non-negative )  eigenvalues,  forms  a  convex 
cone.  We  introduce  a  partial  ordering  with  respect  to 
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the  cone  of  positive  matrices: 

Definition!   A  <  B  means  that  all  eigenvalues  of  B  -  A 
are  non-negative.   Strict  inequality  is  defined  analogously. 

Since  we  have  shown  that  over  a  linear  space  of 
matrices  all  of  whose  elements  have  only  real  eigenvalues 
the  set  of  non-negative  matrices  forms  a  convex  cone, 
it  follows  that  over  such  a  space  this  partial  ordering 
is  transitive,  i.e.  addition  of  inequalities  and  multi- 
plication by  positive  scalars  is  permissible. 

Theorem  3«2l   Let  X  and  Y  be  two  matrices  all  of  whose 
linear  combinations  have  real  eigenvalues  and  suppose  that 
X  is  negative.   Let   jx  be  an  eigenvalue  of  X  +  lY  J 

(3.1)  (X  +  lY)  a  =  iia   . 

Then  the  real  part  of  ^  is  negative. 

Proof!   In  the  symmetric  case  the  proof  follows  by 
taking  the  scalar  product  of  (3.1)  with  the  eigenvector 
a  from  which  it  follows  that    Rep,  =  (a,Xa)  •   In  the 
general  case  consider  the  hyperbolic  equation 

(3.2)  u^  =  Xu^  +  YUy   ; 

it  follows  from  (3»1)  that  the  function 

(3.3)  u'^(x,y.t)  =e<l^'*^*ly>^a   , 

depending  on  the  parameter   "ZT^,      is   a   solution  of    (3»2) 
for   all  values  of    T    ,      Since  by   assvimptlon     X        (X) 
Is  negative,    it   follows  from  formula   (2.i|)    that   the 
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domain  of  dependence  of  the  point  P  =  (0,1)  is  contained 
in  the  half -space  x  <  0  «   In  this  half-space  the  initial 
value  of  our  function  u   defined  in  (3«3)  tends  to  zero 
as  tr   tends  to  plus  infinity,  and  so  do  all  its  deriva- 
tives.  Since  according  to  Theorem  2ol  the  solution  is  a 
continuous  function  of  the  data,   u  (P)  =  e^^   a  must 
also  tend  to  zero  as  "^     tends  to  plus  infinity.  This 
shows  that  the  real  part  of  p,  is  negative© 

Likewise,  if  X  is  non-positive  the  real  parts  of 
the  eigenvalues  of  X  +  iY  are  non -positive* 

Let  X  and  Y  be  a  pair  of  real  matrices  all  whose 
linear  combinations  have  only  real  eigenvalues.  Accord- 
ing to  Theorem  3»1,  the  largest  eigenvalue  X___  of 
X  +  sX  is  a  convex,  the  lowest,   ^t  *   is  a  concave 
function  of  s.   If  in  addition  X  is  positive,  both 
the  smallest  and  largest  eigenvalues  of  Y  +  aX  tend 
to   +  CD  as   s  tends  to   +  oo  »  A  convex  or  concave 
function  with  such  a  behavior  at   +  oo  is  necessarily 
mono tonic.  Thus  we  have  shown  that  the  largest  and 
smallest  eigenvalues  are  mono tonic  matrix  functionso  We 
shall  now  prove  a  more  general 

Theorem  3.3*   Let   ^X^   be  a  linear  space  of  real 
matrices  over  the  reals  such  that  all  eigenvalues  of 
each  X  are  real.   Then  the  eigenvalues  X^(X),^(X), 
,,,,\  (X)   are  mono tonic  functions  in  the  sense  of  the 
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partial  ordering  Introduced  before. 

For  X  and  Y  symmetric,  the  monotonicity  of  the 
eigenvalues  follows  from  the  classical  min-max  principle 
of  Fisher,  [51,  see  e.g.  Riesz-Nagy,  [19],  p.  236,  for 
the  j    eigenvalue; 

Proof:  We  have  to  show  that  the  eigenvalues  X. 
of  Y  +  sX,Y  In  \xj     and  X  positive,  are  mono  tonic 
functions  of  s.  The  eigenvalues  7i(s)  of  Y  +  sX 
are  branches  of  an  algebraic  function  of  s  and  so  have 
only  a  finite  number  of  singularities.  We  shall  show 
that  at  every  regular  real  point  Sq  >  ^  ^^o'  ^^  non- 
negative.  Then  we  can  deduce  that  each  branch  of  X{s) 
is  non-decreasing  in  any  interval  of  regularity.  From 
this  it  follows  that  every  branch  of  X(s)   is  strictly 
increasing;  for  otherwise  X  would  be  constant  in  an 
interval  and  therefore,  being  analytic,  would  be  constant 
for  all  s.   But  since  X  is  positive,  all  eigenvalues 
of  Y  +  sX  tend  to  infinity  with   s  .   Since  9^  is 
continuous  at  the  possible  branch  points,^  this 

ij.  Actually,  7i(Y  +  sX}  is  regular  on  the  whole  real  axis. 
For  once,  the  proof  of  this  result  in  the  symmetric  case 
due  to  Rellich,  [l8],  can  be  carried  out  without  altera- 
tion to  the  more  general  situation.  Let   Sq  be  a 
singular  point  on  the  real  axisj  then  ?i(s)  has  a 
Puiseux  expansion  in  fractional  powers  of  (s-s^)  • 
Since  all  branches  of  7i(s)  are  real  valued  for  s 
real,  there  can  be  no  fractional,  only  integer,  powers 
present  in  the  expansion* 
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completes  the  proof  of  Theorem  3»3,  except  for  showing 
that  \   (Sq)  >  0  for  regular  real   s--,  .   Assume  now  to 
the  contrary  that  at  a  regular  real  point  Sq  --  which 

r 

to  simplify  the  notation,  we  take  to  be  zero  --  \   is 
negative,  say  X  (0)  =  -  a,  a  >  0  »  Then  we  have, 
expanding  about  s  =  0  , 

7l(Y+sX)  =  7i(Y)  -  as  +  0(  ls|^)   , 

Since   s  =  0  is  a  regular  point,  the  above  relation 
holds  for  complex  values  of  s  as  well.   Put   s  =  it 
and  multiply  the  relation  by  i  j   we  get 

?t(-tX+iY)  =  iX(Y)  +  at  +  O(t^)  . 

This  shows  that  for   t  positive  and  small  enough  the 
matrix  -  tX  +  iY  has  an  e i gen value  with  positive 
real  part.  But  since  X  was  assximed  positive,  this 
contradicts  Theorem  3«2« 

We  indicate  now  another  possible  approach  to  the 
results  presented  in  this  section: 

Let   p(4»  V  >f^)     be  a  form  of  degree   n  in 
Cf  *i  »X,  with  the  following  property;   The  coefficient 
of  tJ^     is  one,  and  for  each  fixed  real  choice  for  ^ 
and  ^  ,      p(C,  \  ,X)  tias  only  real  zeros  in  X  •   If 
A  and  B  are  a  pair  of  S3nmnetric  n  by  n  matrices, 
the  form 

-  23  - 


P(?,  \   A)  =  det  1^  +  y|B  -  "Ml 


has  this  property.   Conversely,  I  conjecture  that  every 
i'orra  with  the  above  property  can  be  so  represented.   In 
support  or  this  conjecture  is  the  fact  that  the  number 
of  coefficients  in  p  --  not  counting  that  of  "K^     — 
is  equal  to  the  number  of  elements  of  A  and  B  if 
one  of  these  matrices  is  normalized  by  an  appropriate 
orthogonal  transformation  to  be  diagonal.   If  true, 
this  representation  theorem  could  be  used  to  deduce  the 
convexity  and  monotonicity  of  eigenvalues  over  a  linear 
space  of  matrices  with  all  real  eigenvalues  from  the 
classical  case  of  symmetric  matrices. 

A  similar  conjecture  has  been  made  by  Chandler  Davis 
(personal  communication). 
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ij.«   Finite  Difference  Schemes. 

In  this  section  \-ie  consider  one  level  explicit 
difference  equations  approximating  the  differential 
equation  (2.1),  i.e»  difference  equations  of  the  form 

(U.l)         u(x,t+^t)  =  Zjl  C^'  u  {X+/\,yt)       . 

This  equation  expresses    the  value  of     u     at   the   point 
X     and  at   the   time      t+/\t      as   a  linear  combination  of 
values  of     u     at  time      t      and  at  points     x+/\^.    ♦      The 
displacements     /^.     will  be    taken  to    depend  linearly 
on     ^t    : 

the  vectors     r.     remain  fixed  as     At      tends   to    zero. 
J  — ^ 

Such  a  difference   equation  is    said  to   be    derived 
from,    or  consistent  with,    the   differential   equation 
(2«1)   if   every   smooth   solution  of    (2*1 )   satlsifies 
(i|..l)   approximately,    i.e.    if  the   two    sides   differ  by 
an  amotint     o(At).      Expand     u(t+At )      and     u(x+A  ^») 
into   a  Taylor   series j    the    left   and   right   sides  of 
(i;,l)  can  then  be  written  in   the   form 

ik*-!    )  u  +  A*  u^   +  O(A^) 

and 


i« 


(i;.i  )  SI  cJ  u  +2Z  c-'  u^^  Aj  +  0(A) 
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Substitute     A  * ^^     ^^^     ^t      ^^  given  by   (2«1)   into 

I  .  .  tt 

(i(.,l    }   and     /2i*^i     ^"o^     Zl\i      ^"*°    ^^'^    ^   according  to 

(4.2).      The   two   differential  operators    (i;*!    )   and 

"  2 

(i|.ol    )   have   to   agree   up  to    terms   of   order     /\^      , 

which  leads   to    the  relations 
(4.I4.)  XZ  C^*  r.   =  A 

J 

The  initial  value  u(XsO)  =  <|i(x)  \miquely 
determines  a  solution  of  the  difference  equation  (Ij-d) 
at  all  times  which  are  integer  multiples  of  At  *  The 
scheme  based  on  (/|»1)  is  said  to  be  convergent  if  the 
solutions  ua   of  the  difference  equation  converge 
uniformly  as  /\t  tends  to  zero  to  the  solution  u  of 
the  differential  equation  (2«1)  with  the  same  initial 
value  <1>  ,   for  all  sufficiently  diff erentiable  data  ^    o 

Suppose  the  initial  function  is  an  exponential: 

I       ix£ 

(Kx)  =  e  ^v  ,     v  =  const* 

The  corresponding  solution  of  (I4.0I)  is  also  an 
exponential,  given  by  the  formula 

(i,.5)        u(x,t)  =  ei=^«  (IZ  cJ'e^^J*^}''  V    , 

i   ^^i^ 
where   t  =  N/\t  »  The  matrix  >   C*^  e   *'   is  called 

the  amplification  matrix  or  generating  function,  and  will 
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be  denoted  by  G(^)  . 

Theorem  i;.!:   The  necessary  and  sufficient  condition 
for  the  convergence  of  a  difference  scheme  is  that  the 
eigenvalues  of  the  associated  amplification  matrix  G- ( ^ ) 
should  not  exceed  one  for  real   ^  • 

This  condition  Is  due  to   v  .   Neumann  (see  e.g» 
[i;],  [11]  or  [15]). 

The  necessity  of  this  condition,  as  already  observed 
in  [l5]  within  the  framework  of  Lp   convergence,  follows 
from  the  principle  of  uniform  boundedness  (see  [l]  )•   For 
let   P  be  any  point  not  on  the  initial  hyperplane,  say 
P  =  (0,1)  ;   UA  (P)   is  a  bounded  linear  functional  of 
the  initial  value  4  j 


UA   (P)  =  La   (t    . 


If     ua  (P)      converges   to      u(P)      for  every     c|i      in   the 

N 
space      C        of     N      times  dlfferentiable   functions,    then 

a  fortiori   the   functionals     La  (i      remain  boiinded  for 

every     cj      in      C      ,      and   so    according  to   the   principle 

of  uniform  boundedness,    the   functionals     L      are   uniformly 

bounded. 

Suppose  now  that  v»  Neumann' s  condition  is  violated, 

i.e.  for  some  real  value  of  ^,G(^)   has  an  eigenvalue  "K 

of  absolute  value  greater  than  one  J 

G(.^)v  =  %v  ,  Ul  >  1   . 
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Consider   the   set   of   data 

According   to    (^.5),      ^A^^A    ~  '^A^^^      ^^   equal   to 

As     ^t      tends  to    zero ,    the      C        norm  of     (J  a       grows 
like      (l/A't)^   »     while      |l'A<l'/\l      tends  to   infinity 
exponentially  in     l/At    $      this   shows    that      1?^.|    >  1 
precludes   the  loniform  boiondedness  of   the   functionals 

The    sufficiency  of  the   condition      Ul   5  1     for 
convergence  will  be   shown  by  expanding   the   given  data 
(j)      in   a  Fourier  series.      This   is   a  classical  procedure  J 
since   the   proof  usually   is  carried  out  under  more 
restrictive   conditions,    we   shall   sketch   the   details 
here? 

Leimna  l+.U      Under  the    assumption  of   theorem  i;«l, 
the   powers   of   the   amplification  matrix     Q      satisfy   the 
inequality^  ^'^ 


|lG^(^)ll  <  const.    +   (Nkl) 
Proofs      Write      G      in  the  form 


n-1 


5  By  the  norm  of  a  matrix  we  mean  the  absolute  value 
of  its  largest  component. 
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G   =  YZ.  ^      +  I CMe      ^        -  1)      I 

according   to    (l4.»3)    this  can  be   written  as 

(i;.6)  G  =  I+2ZcMe      J        -1)=I+H      . 


Clearly, 

l|Hl|=0(^)      . 

Let      S     be   a  matrix  such  that 


S'-'-HS  =  h' 


is  triangular.   Such  an  S  exists,  and  it  can  be  so 
chosen  that  the  norms  of  S  and  of  S~   are  bounded 
by  a  quantity  that  depends  only  on  the  order  n  of  the 
matrices  in  question.   Therefore  the  size  of  H   does 
not  exceed  a  constant  multiple  of  |1h||  i 

ik^l)  IIh'II  <  const   llHll     =  0(k|)      . 

Prom   ( i4.«  6  )  we  have  { 

(i|.,6'  )  g'    =  S"-^GS  =  I   +  h'       I 

I  t 

G   also  is  a  triangular  matrix.   Separate   G   into  a 

diagonal  and  subdiagonal  part: 

(Ij-.S)        g'  =  D  +  T   . 

t 

Since   T   is  the  subdiagonal  part  of  H   its  norm  is 
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bounded  by  that  of   H  ;   so  we  have 
{k.7[)  llTll  «  0  (kl)   . 

The  diagonal  elements  of  D  are  the  eigenvalues  of  G   J 
these  are  the  same  as  the  eigenvalues  of  G  which,  by 
the  assumption  of  Theorem  l\.»l,    do  not  exceed  1   in 
absolute  value. 

Raise  both  sides  of  (l+.S)  to  the  Nth  power}  we 
obtain  a  representation  for   (G  )    as  a  sum  of  products 
of  N  factors  D  and  T  •   Since  D  is  diagonal  and 
T   is  subdiagonal,  all  terms  which  include  n  or  more 
factors  of  T   are  zero.   Since  furthermore  the  eigen- 
values of  D  do  not  exceed  one  in  absolute  value  and 
since  the  norm  of  T   is  bounded  by  (ij..7  ),  we  get  the 
following  estimate  for   (G  ) 

1|(g')^1I  <  const.  +  (Nkl)''"^   . 

Since   G^  is  equal  to   S(G*  )^S"*'^  ,   a  similar  estimate 
holds  for  G   ;  this  proves  Lemma  l].,l    • 

Lemma  h.,2t      Under  the  hypothesis  of  Theorem  l;.l, 
the  linear  functionals  La   relating  the  value  of 
solutions  of  (Ij-.l)  at   P  =  (0,1)   to  their  initial 
values   <1»   are  uniformly  bounded  for  all  values  of  A 
in  the   C^  norm  for  M  sufficiently  large o 

Proof  J   Represent   <i>   in  C   as  a  Fourier  transform^  'i 

6  We  assume  that   5   is  zero  outside  of  an  open  set 
containing  the  domain  of  dependence  of   P  . 
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(1+.9)       c^(x)  =  /il/(^)e^''^d?  . 

The  solution  ua   of  (It-.l)  with  initial  value  ^      given 
by  (ij.»9)  can  be  written  in  the  form 

(4.10)  UA  (x,t)  =  /G^(/\t?)^K^)e^^^d?   ,    t  =  N/\t   . 

This  formula  is  obtained  by  superposition  from  formula 
(I|.,5)«   To  prove  its  validity  we  have  to  show  that  the 
integral  on  the  right  in  {1^*10)   converges.   To  this  end 
we  use  the  estimate  of  Lemma  (l+.l)  which  yields  the 
following  majorant  for  the  Integral  (J4..IO); 

/[const  +  (N^tkl)''"^}  U(?)|d^   { 

at      t   =  N/\t   =  1      this   integral    is   equal   to 

/[const   +    kT"^  }    U(^)|dC      . 

But  this  last  Integral  converges  If   \Jf   is  the  Fourier 
inverse  of  a  sufficiently  dlff erentlable  function  cj)  . 
The  order  k  of  the  required  differentiability  depends 
only  on  n  . 

This  reasoning  shows  at  the  same  time  that  the 

value  of   |ua  1   at   P   is  bounded  by  a  constant 

M  I 

multiple  of  the   C   norm  of  (|)  ,   for  M  large  enough* 

This  completes  the  proof  of  Lemma  li.,2. 

Next  we  verify  that  for  an  exponential  initial 

function  (^  =  e   v  ,   the  corresponding  solution  ua 
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of  the  difference  equation  (!|..l)  does  tend  to  the  exact 
solution  u  of  the  differential  equation  (2.1 )•  Accord- 
ing to  (l4.»5),  UA  (x,t)  is  given  by  the  formula 

(il.ll)       UA  (x,t)  =  Q^{/\t^)e^^^   ,     N^t  =  t   . 

Using  the  definition  of  G   and  relation  (14..3),  we  can 
write  G(/\tC)   as 

G(At^)  =  I  +IZ  C^(e   J   -  1) 

=  I  +  i  H  cJ  /\^^   +  O(At^) 

=  I  +  i^t  IZ  C^rj  •  ?  +  O(At^)   . 

Raise  both  sides   to   the     Nth     power,      N  =  t//\    ,      and 
take   the   limit   as     ^  t     tends   to   zero.      Using   the 
classical   rule 

lira        (I   +  ^t  a  +  0(At^))^'^At   =  e'^^*      , 

At-->o 

we   obtain  for    the  limit  on  the   right   in   (1].«1)   the  value 

l[tT~  C-^r.  -5   +  xC] 
e  "^  V     . 

According  to  the  relation  {U.»k.)   between  the  coefficients 
of  the  difference  and  differential  equation,  this  is 
equal  to 
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which  Is,  it  is  easy  to  check,  an  exact  solution  of  the 
differential  equation  (2.1)  taking  on  the  assigned 
initial  values. 

This  shows  that  the  solution  of  the  difference 
equation  with  prescribed  exponential  initial  data  tends 
to  the  corresponding  solution  of  the  differential 
equation.   By  the  principle  of   superposition  the  same 
is  true  for  a  solution  with  finite  linear  combinations 
of  exponentials  as  initial  data.   Since  the  exponential 
functions  span  the  space   C   ,   and  since  furthermore 
the  functionals   La   were  shewn  in  Lemma  1^,2   to  be 
uniformly  bounded,  it  follows  by  a  standard  argument 
(Banach-Steinhaus )  that   L/\(i   converges  to   u(P)   for 


^      in   C^  .   This  cor 
l+.l. 


I       k 

all  (^   in  C   ,   This  completes  the  proof  of  Theorem 


Theorem  I4..I  is  an  exact  analogue  of  Theorem  2.1j 
V.  Neumann's  condition  plays  the  same  role  for  difference 
equations  as  the  condition  of  hyperbolicity  for  differ- 
ential equations. 

If  we  want  to  use  v«  Neumann's  condition  to  verify 
that  some  proposed  difference  scheme  converges,  we  have 
to  check  that  the  eigenvalues  of  the  amplification 
matrix  G(^)   do  not  exceed  one  in  absolute  value  for 
all  real  values  of   ^  .   This  is  not  always  an  easy 
task,  as  any  numerical  analyst  who  has  tried  a  fairly 
complicated  case  can  testify.   Therefore  there  is  room 
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for  more  directly  verifiable  stiff Iclent  conditions  for 
convergence.   Such  a  one  has  been  given  by  Prledrlchs, 
see  [6]  or  [l5J  for  symmetric  difference  schemesj  our 
aim  In  this  section  Is  to  extend  It  to  not  necessarily 
symmetric  schemes. 

Theorem  i|,2j   If  all  linear  combinations  with 
real  coefficients  of  the  coefficients   C"^   of  the 
difference  equation  (L|..l)  have  only  real  eigenvalues, 
and  If  each  C"^   Is  non-negative,  then  the  scheme  based 
on  (i;«l)  Is  convergent. 

We  shall  show  that  under  the  assumptions  of  theorem 
1|.»2  all  matrices  of  the  form 

z.   complex  numbers  of  modulus  one,  have  eigenvalues 
not  exceeding  one  In  modulus.   Since  the  amplification 
matrix  G   Is  of  this  form,  with  z  .  =  e   "^   ,   this 
shows  that   v,  ITeiamann's  condition  Is  satisfied. 

The  symmetric  case  is  easy  to  handle}  let  "k     be 
an   eigenvalue,   v   the  corresponding  normalized  eigen- 
vector; 


^v  =  X_  z  .C-^v   . 

J 

Take   the   scalar  product   of  both  sides  with     v    ; 


■XCvjv)   =  2^  z..(cJv,v)      , 
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The  absolute  value  of  the  left  side  is  \\\    ,      while 
that  of  the  right  does  not  exceed   >   (C'^v,v)  ;   here 
we  have  used  the  fact  that   C'^   is  non-negative.   But 
according  to  (Ij..3)»  this  last  expression  is  equal  to 
(v,v)  =  1  . 

We  turn  now  to  the  nonsyininetric  case.   Consider 
the  matrix 

and  replace   I  by  ">   C"^  ,   as  in  (i|.«3),  and  separate 
real  and  imaginary  parts.  We  get 

where  x.   and  y.   are  the  real  and  imaginary  parts 
of  z.  .   Since  x.   cannot  exceed  1  ,   the  real  part 
of  the  above  matrix  is  a  linear  combination  of  non- 
negative  matrices  with  nonpositive  coefficients  and 
is  therefore  by  our  convexity  theorem,  itself  non- 
positive.   Then,  according  to  theorem  (3«2),  the  real 
part  of  the  eigenvalues  of   >   z.C*^  -  I   are  also  non- 
positive.   So  we  have  shown J 

The  real  parts  of  the  eigenvalues  of   "^   z  .C "^ 

J 

do  not  exceed  one»   But  since  the  eigenvalues  of 
y       z.C^    are  homogeneous  functions  of  degree  one, 
this  shows  that  the  moduli  of  the  eigenvalues  cannot 
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exceed  1  either* 

It  is  easy  to  show  that  the  conditions  in  theorem 
i|«2  on  C"^   are  not  only  siofficient  but  necessary  for 
y        z.C''   to  have  all  its  eigenvalues  in  the  unit  circle 
for  all  possible  choices  of  z.  ,  |z^|  =  1  «   In  case 
that  the  vectors   r.  do  not  lie  on  a  lattice,  it  can 
be  shown  that  v*  Neumann' s  condition  implies  that 
>   z.C**  has  no  eigenvalues  of  modulus  greater  than 
one  for  all  possible  choice  of  ^a    »    \zA    =  1.    (see  [9] 
for  a  proof )e   Therefore  in  this  case  the  conditions 
contained  in  theorem  i|.. 2  are  necessary  as  well  as  suffi- 
cient for  the  convergence  of  a  difference  scheme. 

In  most  difference  schemes  the  coefficients   G*' 
are  linear  combinations  of  the  coefficients  of  the  hyper- 
bolic differential  equation  which  the  scheme  is  supposed 
to  approximate.   In  such  a  situation  the  first  assumption 
of  Theorem  l\.c2   --  about  the  reality  of  the  eigenvalues  — 
is  automatically  fulfilled,  since  the  coefficients  of 
the  given  hyperbolic  equation  have  this  property. 

In  the  symmetric  case  Prledrichs  has  proved  more; 
he  has  shown  that  \xader  the  hypothesis  of  Theorem  Ij.*2 
the  difference  operator  (I4..I)  relating  the  value  of  u 
at   t  to  its  value  at   t  +  ^t   does  not  increase  the 
Lp  norm  with  respect  to  the  space  variables  I 

/u^(x,t+/\t)dx  <  /u^(x,t)dx   • 
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Furthermore  he  has  shown  that  even  If  the  coefficients 
C*^   depend  on  the  variables  x   and   t  ,   we  have 

ru^(x,t+^t)dx  <  (l-K)(/\t))  /  u^(x,t)dx   , 

provided  that   C*'   are  Llpschitz  continuous.   This 
guarantees,  it  is  easy  to  show,  the  convergence  of  the 
difference  scheme. 

We  shall  now  present  some  applications  of  Theorem 
1+.2. 

i)  Convert  the  differential  equation  (2,1)  into  a 
difference  equation  by  replacing  the   x   derivative  by 
a  symmetric  difference  quotient! 

u(x+5j^ej^)  -  u(x-5j^e^) 

V ^^ 

where   e.  is  the  unit  vector  in  the  x*^   direction,  and 

J 

by  replacing  the  time  derivative  by  a  forward  difference 
quotient  J 

u(x,t  +  ^t )  -  u(x,t) 

U  .  ^r^ , 

where  u  is  the  mean  value  of  u  at  the   2m  points 
x  +  6^e^      t 

^[u(x  +  \ej^,t)  +  u(x-5j^e^,t)  j 
u(x,t)  =^^ ^ 

The  resulting  difference  scheme  is  of  the  form  (ij.«l), 
with  the   2m  displacement  vectors   +  5,  e,  ,  and  with 
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coefficients 

+k        2in       —  2o, 

Applying  Theorem  ^,2   shows  that  this  difference  scheme 
is  stable  if  /\t      is  chosen  so  small  that  &-^^m/\t 
exceeds  the  absolutely  largest  eigenvalue  of  A   •   This 
is  am  improvement  of  a  result  due  to  Fritz  John  who  has 
shown  (in  an  unpublished  note)  that  the  above  difference 
scheme  is  stable  if  /\^t   is  chosen  to  be  small  enough, 

ii)  Susan  Hahn,  in  [9],  has  investigated  slmplicial 
difference  schemes,  i.e.,  schemes  in  which  the  number 
of  displacement  vectors  r.   is  one  greater  than  the 
number  of  space  variables  x   and  In  which  the  displacement 
vectors  are  the  vertices  of  a  simplex.   She  has  shown 
that  if  such  a  scheme  satisfies  the  classical  necessary 
condition  for  convergence  due  to  Courant,  Priedrlchs  and 
Lewy,  [2],  i.e.,  if  the  domain  of  dependence  of  the  point 
P  =  (0,1)   lies  in  the  simplex  formed  by  the  displacement 
vectors  r.  ,   then  the  coefficients   C.   are  non-negative. 
Thus  she  was  able  to  conclude  —  using  Priedrlchs'  theorem 
in  the  symmetric  case  and  Theorem  i;«2  in  the  general  case  — 
that  for  a  simpllcial  scheme  the  Courant -Priedrlchs -Lewy 
condition  is  sufficient  as  well  as  necessary  for 
convergence. 


-  38  - 


BIBLIOGRAPHY 


[1]      Banach,    S. ,    Operations  Llnealres,    1932.   Warsaw. 

[2]      Courant,   R. ,    Frledrlchs,   K.    C,    and  Lewy,    H,, 
tJber   die   partlelle   Dlff erenzlalglelchurgen  der 
mathematischen  Physlk,    Math.    Annalen,    Vol.   100, 
1920,    pp.    32-74o 

[3]      Courant,    R.    and  Lax,    A.,   Remarks   on  Cauchy's 

Problem,   Comm.    on  Pure   and  Applied  Mathematics, 
Vol.    5,    1955,    pp.    i;97-502. 

[I|.]      Eddy,   R.    P.,    Stability  in  the  nujierical    solution 
of  initial   value   problems    In  partial   differential 
equations.  Naval   Ordnance   Lab,    Memo.    10232. 

[5]   Fisher,  E. ,  Uber  quadratische  Formen  mit  reellen 
Koefficienten,  Monatshefte  fur  Math,  und  Phys., 
Vol.  16  (1903),  pp.  23i;-2i+9. 

[6]  Friedrichs,  K.  0.,  Symmetric  Hyperbolic  Linear 
Differential  Equations,  Comm.  Pure  and  Applied 
Mathematics,  Vol.  7,  195k,    PP»  3i|5-392. 

[7]   Garding,  L.,  Linear  hyperbolic  partial  differential 
equations  with  constant  coefficients.  Acta  Math., 
Vol.  b3,  1951,  pp.  2-62,    ~ 

[8]   Hadamard,  J.,  Le  probl^e  de  Cauchy  et  les 
equations  aux  deriv^es  partlelles  lingairis 
hyperboliques  (Herman,  1932). 

[9]   Hahn,  S. ,  Dissertation,  New  York  University,  1957« 

[10]   Holmgren,  E. ,  Uber  Systeme  von  llnearen  partiellen 
Differential  glelchurgen,  'Ofversigt  Konglo 
Vetens.-Akad.  Forho ,  Vol.  58,  1901,  pp.  91-103» 

[11]   Hyman,  M.  A.,  O'Brien,  G.  G. ,  Kaplan,  S. ,  A  Study 
of  the  niunerical  solutions  of  partial  differential 
equations.  Journal  of  Mathematical  Physics,  Vol.  29, 
1951,  pp.  223-251. 

[12]   John,  F.,  On  linear  partial  differential  equations 
with  constant  coef f icients^  Comm.  on  Pure  and 
Applied  Mathematics,  Vol.  2,  191+9,  pp.  209-253» 

[13]   John,  Fritz,  Plane  Waves  and  Spherical  Means, 
New  York,  Intersclence,  1956. 

-  39  - 


[li;]   LsLX,  A.,  On  Cauchy' 3  problem  for  partial  differ- 
ential equation  with  multiple  characteristics, 
Commo  Pure  and  Appl.  Math.,  Vol.  9,  1956, 
pp.  135-168. 

[l5]   Lax,  P.  D. ,  Rlchtmyer,  R.  D. ,  Survey  of  the 

stability  of  linear  finite  difference  equations, 
Comm.  Pure  and  Appl.  Math.,  Vol.  9,  195o» 
pp.  267-293. 

[16]   Lax,  P.  D. ,  Asymptotic  solutions  of  oscillatory 
initial  value  problems,  Duke  Math.  Journal,  to 
appear  in  December  19^7» 

II?]   Plancherel,  M. ,  Polya,  G.,  Fonctions  entieres 
et  integrales  de  Fourier  multiple.  Comment arii 
Matematici  Helretici,  Vol.  9,  1936-37,  pp.  22i|-2i4.8. 

[l8]   Rellich,  F.,  Storun^stheorie  der  Spectralzerlegung, 
Math.  Annelen,  Vol.  113  (1936),  pp.  600-619. 

[19]   Riesz,  F. ,  Nagy,  Szo-B. ,  Lejons  d' analyse 

fonctionnelle,  Akadlmiai  Kiadd,  Budapest,  1952. 


-  ko  - 


y 


Date  Due 

JAM  30 -es 

APR:i^ 

IS17 

■       K:  r-5 

— *^  "^ 

.::'^ 

:.:a,S24^ 

~^^  ^— 

^ 

^'-•s^ 

,-.u3  2  7 

p- 

!          1 

iEPi 

1  1330 

mzrM 

'  - 

-aA 

UAR   9 

'au 

tSP  e;  • 

-i" 

DEC  ITl-e 

1 

1 

JAN  S 

;« 

j)/. 

•« 

5 

""jT-T     * 

•ftY 

NOV  U  ■ 

•t 

Mir    '      1 

MAR  6   -'9^ 

mt  2  2  ■ 

*f*-r     r»  9  10' 

C 

OCl     ^  c 

m 

PHINTeO 

IN  U.  S.  A. 

/  rJYU 

NYO 
797l| 


Lax 

Differential   equations. 
Difference   equations   and 
ma t r ix   theory. 


